The renormalized jellium model for spherical and cylindrical colloids 
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Starting from a mean-field description for a dispersion of highly charged spherical or (parallel) rod- 
like colloids, we introduce the simplification of a homogeneous background to include the contribution 
of other polyions to the static field created by a tagged polyion. The charge of this background is self- 
consistently renormalized to coincide with the polyion effective charge, the latter quantity thereby 
exhibiting a non-trivial density dependence, which directly enters into the equation of state through 
a simple analytical expression. The good agreement observed between the pressure calculated using 
the renormalized jellium and Monte Carlo simulations confirms the relevance of the renormalized 
jellium model for theoretical and experimental purposes and provides an alternative to the Poisson- 
Boltzmann cell model since it is free of some of the intrinsic limitations of this approach. 



I. INTRODUCTION 



In the realm of soft matter, a quantitative description 
of systems with a non vanishing density of mesoscopic 
constituents (colloids) is a difficult task whenever long 
range Coulomb interactions are present It is cus- 
tomary to introduce a suitably defined Wigncr-Scitz-likc 
cell to render the situation tractable. This considerable 
simplification allows for the computation of various ther- 
modynamic quantities (see e.g. 0, [|| and [J, Q for more 
recent accounts). Transport properties may also be de- 
rived [1, 0| ■ Initially borrowed from solid state physics, 
the concept of a cell nevertheless appears fruitful to de- 
scribe the phase behaviour of liquid phases (see e.g. Q). 
When a more microscopic information such as effective 
interaction is sought, there is however to date no evidence 
that the cell picture provides accurate answers, through 
the approximate procedures that have been proposed to 
infer solvent and micro-ions averaged inter-colloid poten- 
tials of interaction In the present paper, we adopt 
a more "liquid-state" viewpoint Q to describe the lo- 
cal and global screening properties of microscopic ions 
around highly charged rod-like or spherical colloids, tak- 
ing due account of the finite density of colloids. Our ap- 
proach bears a strong resemblance with a Jellium model 
put forward by Beresford-Smith, Chan and Mitchell [l(|, 
with the important difference that the jellium under con- 
sideration here is the "renormalized" counterpart of that 
studied in (lQ [. A preliminary account with emphasis 
on the sedimentation of charged colloids has been pub- 
lished in 11[ and we note that our approach has been re- 



cently tested with some success on liposome dispersions 
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The article is organized as follows. The model is de- 
fined in section [II] and illustrated in section IIIII where 
salt-free suspensions are considered and the numerical 
procedure exemplified with charged spheres. The cylin- 
drical geometry is also addressed which allows to discuss 
the fate of classical Manning-Oosawa condensation phe- 
nomenon [lj| within the present framework. The ef- 
fects of an added electrolyte are investigated in section HVl 
and concluding remarks are presented in section [Vj As 
will become clear below, our framework provides a pro- 
cedure to incorporate self-consistently charge renormal- 
ization into the classical Donnan equilibrium description 
of suspensions. 



II. THE MODEL 

When immersed in a polar solvent, mesoscopic "par- 
ticles" release small ions in the solution that, together 
with other micro-ions resulting from the dissociation of 
an added salt, form an in- homogeneous cloud around 
each colloid. Within the mean-field approximation to 
which wc restrict here, neglect of micro-ionic correlations 
allows to relate the local density of micro-species of va- 
lency z.i at point r to the local electrostatic potential 
ip(r) through rii(r) oc exp(— Zie(3Lp(r)), where e is the 
elementary charge and = kT is the thermal energy 
P, [l?]]- For any position of the N colloids present, one 
needs to solve the resulting Poisson-Boltzmann equation 
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from which the electric potential follows. This potential 
may then be inserted into the stress tensor [l[ to com- 
pute the force acting on the colloids. Such a procedure, 
which makes explicit use of the separation of time scales 
between colloids and micro-ions, opens the way towards 
a complete description of the statics and dynamics of 
the system, with e.g Monte Carlo or Molecular Dynam- 
ics techniques to treat the colloidal degrees of freedom 
[3 . This treatment is however numerically particu- 
larly demanding and much insight is gained from further 
approximations that map the original problem with TV 
colloids onto a one-colloid situation. The cell model is an 
option, where the finite density of colloids is accounted 
for by an exclusion region. Wc propose here an alterna- 
tive, free of some of the limitations of the cell model, that 
is equally simple to implement and solve. 

A given colloid with bare charge Z^^e, assumed pos- 
itive, is tagged and fixed at the origin. The suspend- 
ing medium (solvent treated as a dielectric continuum 
with permittivity e) is taken as infinite, with a mean 
colloidal density p. Following [nj, the colloids around 
the tagged particle are assumed to form an homogeneous 
background, with charge density Z^^ep, so that the 
electrostatic potential around the tagged colloid fulfills 
Poisson's equation 



V 2 (y5 = 



-fieznp 



(i) 



where the summation runs over all micro-species and 
the concentrations n° are determined either from electro- 
neutrality in the no salt case or from the osmotic equilib- 
rium with a salt reservoir in the semi-grand-canonical sit- 
uation, both addressed here. At large distances (r — > oo), 
the term in brackets on the rhs of ([T]) vanishes which im- 
poses a value (foe for the potential far from the colloid, 
that may be called a Donnan potential. The key point 
in our approach is that unlike in [To| . Zback ^ Zb&re'- the 
background charge is not known a priori but is deter- 
mined self consistently as explained below. 

To illustrate the methodology, we consider a spherical 
colloid of radius a. When r — > oo, we may linearize 
(fTj) around ifoo, which results in a Hclmholtz equation 
indicating that 
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where the characteristic distance « -1 is given by 
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and £b = /3e 2 /e is the Bjerrum length. 

For very low bare charges, the solution ([2]) holds for 
all distances with Z e g = Zbare, and one can consider that 
■Zback = ■Zbarc- However, typical colloidal charges are 
such that Zbare^s/a 3> 1, a regime for which countcri- 



ons become strongly associated with the colloid and the 
charge renormalization effects [2(| [U [22| can not be ig- 
nored. The counterion condensation strongly affects the 
electrostatic far-field so that the large distance signature 
involves an effective charge [Z c s in Eq. ([2])] which signifi- 
cantly differs from the bare one. As a result of non-linear 
screening, one has Z e g <C ^barc whenever Z^^Ib/ci 1. 

At this point, the effective charge arising in @ is a 
function of both the background and the bare charge, 
other parameters being fixed: Z c g = Z c ff(Zb ac k, ^baro)- 
As far as a static description is pursued, for sufficiently 
strongly charged colloids the bare charge is an irrelevant 
quantity far enough from the tagged colloid, and we de- 
mand that Zback coincides with Z e g, which best char- 
acterizes the background charge resulting from smearing 
out the other colloids contribution. We therefore enforce 
the self-consistency constraint 
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(4) 



to compute the a priori unknown background charge. As 
we shall see below, this condition is readily implemented 
numerically and for a given ^bare, leads to a unique value 
for Zback = Z e ff. This value is density dependent, which 
is also the case of the inverse screening length n. Indeed, 
ifioo depends on p (23j through the electro-neutrality con- 
dition 



-^back p e 



;e exp (-^eziipoc) = 



(5) 



which translates into a p dependence for k. Considering 
now two colloids in the weak overlap approximation (i.e. 
not too close), the effective potential of interaction will 
take a DLVO form [l], [l7| with effective parameters k and 



The procedure outlined here incorporates non linear 
screening together with finite p effects. It is best suited 
to describe low density systems since the colloid-colloid 
pair distribution function g cc is implicitly considered to 
be unity for all distances. This reduction, which has non 
trivial consequences, is certainly of little relevance for 
high density suspensions for which the cell model is pre- 
sumably a better approximation. 

Before illustrating the method, we briefly consider the 
pressure in the system, that is given by the densities of 
micro-ions far from the tagged colloid: 

f3P = ^2 »i(oo) = n° exp (-j3ezicpoo) ■ (6) 

i i 

The colloidal contribution is explicitly discarded [ill ]. 
This is well justified in the low salt limit, which is a 
regime of counterions dominance provided that Zb are 3> 
1, which is easily achieved in practice. 
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III. THE NO SALT LIMIT 

A. Spherical colloids 

The simplest situation to investigate is that of de- 
ionized suspensions (no salt). For simplicity, we consider 
counterions as monovalent. From (fl]) it follows that the 
dimensionless potential <f> = f3etp obeys the equation 



dr 2 



2d<j> 
f dr 



„ -Zback^B , a 

6-q (e* 



1 



(7) 



where a is again the radius of the tagged particle from 
which the dimensionless distance r = r/a is defined, and 
i] = Airpa 3 /3 is the volume fraction. The boundary con- 
ditions are 



0^0 



Zi 



bare 



dr 



for f 
for f 



1. 



(8) 
(9) 



In writing Q , use has been made of the (global) electro- 
neutrality constraint n°_e cxp((/) oc ) = Z^^pe with the 
choice (f>oo = 0. For all values of Zback, the far-field of <f> 
is governed by k such that 



(rea) 2 = 3i] 



Zback^B 



(10) 



The above system is solved following similar lines as 



24j . We summarize here the main steps. In prac- 



tice, equation ([7]) is solved numerically for a finite sys- 
tem f G [1, -R], where i? needs to be large enough (that is 
naR 3> 1 but note that k is not known initially but fol- 
lows once the background charge is known), a) The first 
and important step is to rephrase the boundary value 
problem at hand as an initial value problem with bound- 
ary conditions 4>'{R) = (to ensure electro-neutrality) 
and varying (f>(R). The volume fraction r/ is fixed and 
the background charge Zback first assigned a guess value, 
to be modified itcratively (see below). If (f>(R) is small 
enough, the system then admits a solution, b) From this 
solution, one computes <p'(r = 1) to know the corre- 
sponding bare charge, c) Changing <ft(R) [2{|, the tar- 
geted value 4>'(f = 1) = — Zbarc^s/a is readily found by 
iteration, d) The screening quantity k is subsequently 
computed from (flU|) and the effective charge associated 



to the particular couple (Zback, Zb are ) is deduced from 
the large f behavior of <fi [e.g one needs to observe a 
clear-cut plateau for [(j)(r) — 4>(R)]e K,ar f plotted as a func- 
tion of f in the range 1 <C f < R\. The first itera- 
tion ends here, and the procedure is repeated with the 
Z e g obtained as the next trial value for Zback- Alterna- 
tively, one may sample several trial values for Zback an d 
plot Z c g versus Zb ac k- As ma y be observed in Fig. [T] 
where such a plot is displayed, the dependence of Z e g on 
Zbare is very mild, which means that convergence towards 
-Zback = Z c g is achieved in a few steps. In the (artificial) 



limit where Zback — * 0, the problem at hand reduces to 
an unscreened one (governed by Laplace equation) with 
solution <fr = Z^qIb/t: there is therefore no renormal- 
ization of effective charge so that Z c r — > Zbare (see Fig. 
[IJ. The inset shows how the self-consistent background 
charge is determined, the other points being unphysical. 

In the limit of a diverging bare charge, the procedure 
is well behaved and yields a finite self-consistent effec- 
tive charge. From the previous discussion, we expect Z e g 
to diverge at small Zback, which is indeed the case (not 
shown) . 

Once the physical solution to the problem has been 
located (inset of Fig. [1]), various quantities such as the 
pressure may be computed. In the remainder, we will 
use the terms "effective" and "background" charges to re- 
fer to the self-consistent solution as obtained in Fig. Q] 
Z c ff = Zback is therefore a function of Z^ alc and vol- 
ume fraction (possibly also salt concentration, see sec- 
tion [IV]). For a particular density, this function is shown 
in Fig. O After the initial linear regime, where no 
renormalization takes place, the effective charge slowly 
reaches a saturation plateau as .Zbarc - * oo. For this spe- 
cific density (77 = 10 -2 ) the effective charge saturates at 
ZcS^b/o- — 6.6. The saturation phenomenon observed 
here is strongly reminiscent of that observed in the clas- 
sical Poisson-Boltzmann approach (either in a cell, or in 
an infinite medium pH I22I |26| ) . To assess quantitatively 
the possible difference with the results obtained within 
the cell model, we compare in Fig. [3] Z c g derived in 
the cell 0, [24[ to its renormalized jellium counterpart. 
Both charges differ by a notable amount for r] > 10~ 3 
while the agreement at very low density is meaningless, 
and follows from a divergence of the saturation effective 
charge in both models: non-linear effects disappear when 
77 — > 0, so that Z c ff — > Zbarc- This is a peculiarity of sys- 
tems with colloidal spheres and counterions only, and it 
turns out that the behavior of charged cylinders is quite 
different, see below. 

Under the de-ionized conditions studied here, the pres- 
sure takes the simple form (3P = Z c gp, whereas the cor- 
responding expression in the cell model is less explicit 
and does not directly involve the effective charge. Re- 
markably, although there is a significant difference be- 
tween the effective charges calculated within the Poisson- 
Boltzmann cell model (PBC) and the renormalized jel- 
lium, the pressures calculated using the two models are 
identical for rj < 0.1, Fig. [5] A similar agreement is 
found at saturation [ll|. We add here that a compari- 
son between the renormalized jellium equation of state 
and the "exact" results of the primitive model obtained 
using the Monte Carlo computations has been reported 
in [ll| , with excellent agreement (the corresponding den- 
sity range is quite low such that PBC and renormalized 
jellium predictions agree). 

Before concluding this section, we emphasize that one 
must carefully check that the results obtained do not de- 
pend on the particular value chosen for the cutoff i?, e.g. 
by repeating the analysis with an increased cutoff. 
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B. Cylindrical colloids 



Consider a nematic phase of parallel infinite rods (L — > 
oo ) with line charge Abarc (therefore no positional order 
in the plane perpendicular to the main axis) . We may re- 
peat the previous approach, tagging a given rod of radius 
a and model the effects of the other rods by an homo- 
geneous background, with line charge Aback- While the 
definition of k is unaffected compared to the spherical 
case, the far-field potential now takes the form: 



r = 



2A 



>cff«B 



Ko(Kaf) 
KaKi(na) ' 



(11) 



where K (resp. K{) denotes the zeroth (resp. first) 
order modified Besscl function of the second kind. In 
the spirit of the consistency requirement of section [Til 
we impose Aback = A c ff where again, A e ff follows from 
the large distance behavior of the solution to Poisson's 
equation with background charge Aback : 



d 2 <f> ld</> 

-jzj + = 4r?A b ack*B 

ar z r ar 



3* - 11 



(12) 



Here the volume fraction is r\ = ira 2 n 2D where n 2D is the 
mean surface density of rods (in the plane perpendicular 
to their axis). The boundary conditions are the same as 
© and @, and the numerical method identical to that 
used in the spherical case. 

The effective charges calculated using the cell and the 
renormalized jcllium models are compared in Fig. [5] for 
Afcare-^B = 1. The inset corresponds to the saturation 
regime where Abarc is very large (Abarc — ► oo). We ob- 
serve a substantial disagreement between the two effec- 
tive charges. On the other hand, in the small bare charge 
regime where non-linear effects are not at work, both 
quantities coincide (not shown), which is a signature of 
whenever non-linear effects come into play (i.e. outside 
the small bare charge linear regime, which is the case 
for both figures). Be yon d these differences, Manning- 
Oosawa condensation [l5l . Il6j , which is a key feature of 
2D electrostatics, is shared by both PBC and renormal- 
ized jellium models. As the colloid density is decreased 
(r) — > + ), the effective charge becomes independent of 
the bare one, provided Abarc exceeds the critical threshold 
1/Ib- This feature is illustrated in Fig. [Sj At the satura- 
tion plateau and again for 77 — > + , one has A er r^B — 0.47, 
a value that will be refined below. 

To be more quantitative, it is furthermore natural to 
compare the corresponding functional forms of effective 
charges versus bare ones, and versus density in both PBC 
(where it can be computed analytically) and renormal- 
ized jellium models, where this information is accessed 
numerically To this end, we reconsider the analytical 
results obtained in [22|, [24[ where the effective charge in 
the cell model following Alexander et al. prescription Q 



reads: 

Acfi^B 



-Kl B aR ws {Ii(Kp B RwB)Ki(KpBa) 
- hiKpB^K^KpBRw.)}, (13) 



with standard notation for the Bcssel functions. Here 
R ws = ar/ -1 / 2 is the radius of the cell and Kpb is the 
inverse screening length related to the micro-ionic den- 
sity at the cell boundary [22| , which can be computed 
explicitly from the analytical solution of (2j. After some 
algebra, we find, to leading order in density that when 
Abarc > 1/£b 



V2 



h(V2) 



,/ (V2) + V2/ 1 (y / 2) + / 2 (^) 
(2£-logfo)) s 



(14) 

where £ = Ab a re/(Abare — 1/^b). We note that the lead- 
ing term v / 2/i(v / 2)/2 ~ 0.63 differs from the value found 
in the renormalized jcllium (~ 0.47, see Fig. [6]). More- 
over, Eq. I|14p also suggests a fitting form to describe 
the saturation plateau in the low density regime of the 
renormalized jellium model: 



Ks£b 



.4 



D 



(C- log 77)2' 



(15) 



The values of A, B and C can be obtained using a numer- 
ical fit. We find that in the saturation limit A ~ 0.471, 
B ~ 16.87 and C ~ 0.843 give an excellent agreement 
with the numerical data. We have also checked that an 
equally good agreement is found at lower bare charges, 
such as Abare^s = 4, but with different values of A, B and 
C. We conclude here that both models are described by 
the same limiting law for low densities, at least beyond 
the condensation threshold. 

It is of interest to resolve the condensate structure once 
the counterion condensation has set in. A useful mea- 
sure of the condensate thickness is provided by the so- 
called Manning radius Rm [13 that has been recently 
worked out in the infinite dilution limit and for low salt 
content 0, H^|: in practice, the integrated charge per 
unit length q(r) around a rod has an inflection point at 
r = Rm, when plotted as a function of logr. This is 
exactly the point where g(i?Af)^s/e = 1. We expect a 
similar behavior for the renormalized jellium, given that 
in the vicinity of highly charged rods, the (largely dom- 
inant) counterion distribution should not be sensitive to 
the difference between a uniform background as in the 
renormalized jellium model, and coions as in the situa- 
tion worked out in [l6[. The lower inset of Figure [7] shows 
that this is indeed the case. In addition, from the analyt- 
ical expressions derived in [l6[ and the fact that the rel- 
evant screening parameter reads here (ko) 2 = 4r)\ c ff£ Bl 
we expect the scaling kRm oc (ko) 1 / 2 , more precisely 



D 1^° -1/4 

Rm oc arj ' exp 



1 



2(Abarc^3 - 1) 



(16) 
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The dependence of Rm on both density and bare charge 
embodied in Eq. (|16[) is fully supported by the numerical 
data, see Fig. [71 

Finally, and much like for spherical colloids, there is 
a good agreement between the osmotic pressure calcu- 
lated using the cell model and the rcnormalizcd jellium 
approximation, in spite of the different effective charges, 
see Fig. El Discrepancies are observed only for volume 
fractions 77 > 0.06 and the agreement seems to be better 
at high charges. 



IV. EFFECTS OF ADDED SALT 



In this section, we consider systems dialyzed against an 
electrolyte reservoir with the monovalent salt concentra- 
tion c s . The corresponding screening parameter is K 2 es = 
8it£bc s . It is convenient to choose the reference potential 
so that micro-ionic densities are n±(r) = c s exp[=F</>(r)], 
where the counterions are assumed to be monovalent. Us- 
ing Eq. ([S]), the potential at infinity becomes 



. , / -^backp\ 

arcsmh — 

V 2c s J 



(17) 



It is important to keep in mind that n±(r) are not the 
physical microion densities, but are the effective (renor- 
malized) quantities satisfing 

J dr[n + (r) - n_(r) + Z ctf p] = -Z haie . (18) 

Since the renormalization does not affect coions, their 
concentration inside the jellium with one colloid fixed at 
r = is 



c + = y I 



(19) 



where it is understood that V denotes the measure of 
a large volume centered at r = 0. The concentration 
of counterions, C_, then follows from the overall charge 
neutrality inside suspension, C_ = C+ + Z^^p. 

Far from colloid, n+(r) saturates at the bulk value n + , 
so that in the thermodynamic limit (V — > 00) 



C+ = n+. 

Similarly, for V — ► 00 

y J dr[n + (r) - n_(r) + Z cS p] = 
which means that 



V 



(20) 



0, (21) 



C + - i J dm-(r) + Z cS p = 



(22) 



The charge neutrality allows us to rewrite Eq. (|2"2")l as 

1 J drn„ (r) = C_ - (Z haie - Z cS )p . (23) 

Eq. (j2"3")l provides a suggestive interpretation of n_ (r) as 
the local density of free (uncondensed) counterions. Far 
from colloid, rt_(r) saturates at its bulk value n_, and 
in thermodynamic limit Eq. (|23p reduces to 



C- = n_ + (Z haIe - Z cS )p. 



(24) 



Eqs. (|20|) and ([24)1 allow us to calculate the ionic con- 
tent inside a suspension dializcd against a salt reservoir. 
This is particularly useful when comparing the results of 
the renormalized jellium model, which is grand canonical 
in electrolyte, with the Monte Carlo simulations, which 
are usually performed in a canonical ensemble. Knowl- 
edge of the assymptotic potential allows us to obtain the 
concentrations of coions and free counterions inside the 
suspension, 



n± 



c s exp(=F0oo)- 



(25) 



These are precisely the densities that govern screening 
within the renormalized jellium, k 2 = At:Ib(p+ +«_). 



A. Spherical colloids 

In spherical geometry, Poisson equation ^ now takes 
the form 



df 2 



2d(j) 
f df 



(K les a) 2 sinh0 — 377 



Z h 



ack^B 



(26) 



We again solve it numerically as a boundary value prob- 
lem in a (large enough) finite cell with vanishing (f>' at 
the boundary, increasing gradually the boundary poten- 
tial from the value 



arcsinh 



3ry^back^s/a 



(Atresa)^ 



(27) 



which corresponds to a vanishing bare charge. 

Linearizing Eq. (|26p around tp^ , it can be seen that at 
large distances the potential takes the form of Eq. (J2J, 
with a screening constant n given by 



(«r 



/ 3r/Z ba ck^B 

V a 



(28) 



For highly charged colloids and typical salt conditions, 
the corresponding density dependence is shown in Fig. 
IH1 while the effective charge (deduced from the condition 
Z c g = Zback) is displayed in Fig. [TUJ When 77 — * 0, both 
quantities coincide with the infinite dilution limit of the 
traditional PB theory, as they should. The increase of k 
with the density of colloids reflects the increasing impor- 
tance of counterion screening. The effective charge shows 
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a non-monotonous behaviour with respect to density. 

To compute the osmotic pressure, we subtract the 
reservoir pressure (2c s fcT) from the expression ((6|). More- 
over, it should be remembered that such a relation only 
provides the ionic contribution to the pressure. In the 
presence of salt and at low colloidal density this con- 
tribution becomes smaller than the colloidal one. The 
vanishing of the microion contribution to pressure is ex- 
ponential in the cell model, while it is algebraic for the 
jcllium. Both models should then strongly disagree in 
the low density limit. To mimic the colloidal contribu- 
tion, we add the ideal gas term pkT to ([6]), so that the 
resulting osmotic pressure reads 

/3n = p+ v /z c V + 4c2-2c a . (29) 

In the no salt case, addition of the ideal term is irrele- 
vant since it is always much smaller than the micro-ionic 
one, provided that Z e g is large enough (this is the case 
for highly or even weakly charged colloids since a £b ) ■ 
Moreover, addition of the ideal gas term breaks the scal- 
ing form valid in the no salt case where a 2 £s/3P only 
depends on j] and reduced charge Z^^b / a- We there- 
fore show the osmotic pressure in Fig. [11] for two values 
of colloid radius, within both the PBC and the renormal- 
ized jcllium frameworks. Apart from the expected devi- 
ations at small densities, one observes compatible values 
at higher rj. 

There exist relatively little simulational data for the 
primitive model with salt, where the bare Coulomb inter- 
actions between all charged species -colloids and micro- 
ions- are taken into account (with still an implicit sol- 
vent}. A reference equation of state with salt is provided 
in [32j, with the simplification of a Wigner-Seitz cell, but 
explicit micro-ions. The simulations were performed in 
canonical ensemble with fixed salt content. The amount 
of added salt is characterized by a ratio of the over- 
all added cation charge to the overall macroion charge, 
0l = C+/(Zbarc p)- We compute the densities C± corre- 
sponding to a given salt content as discussed in section 
ITVl In Fig. [12] the osmotic pressure (3P/ p t is plotted as 
a function of /3l, where pt is the total density of ionic 
species. As in the case of salt-free suspensions the pres- 
sures calculated using the PBC and the renormalized jel- 
lium are in good agreement. 

B. Rod-like colloids 

For completeness, we briefly report here results for 
cylindrical geometry. Unlike salt-free case where A e ff is a 
monotonic function of density, a minimum appears in the 
renormalized jellium curve shown in Fig. 1131 The agree- 
ment between PBC and renormalized jellium at low 77 sig- 
nals the region where the system is salt dominated (the 
colloid density is too low and, consequently, countcrions 
do not participate in screening). Conversely, the inset in- 
dicates the density range where counterions do dominate 



: for X] > 10 -1 , the results become independent of the 
reservoir ionic strength and coincide with those obtained 
in the no salt limit . 

Finally, the pressure ([29]) for cylindrical colloids is 
given by 



4^ B a 2 /?n = 4ry-^ + v /(4A e// £ B r/) 2 + ( Kres a) 4 -( Kl . es a) 2 . 

(30) 

Note that for infinite polyions [L — > 00), the first term 
on the right hand side of Eq. ([30]) vanishes. In Fig. [14] 
we plot the equation of state for polyions of Abare^s = 2. 
One should note a strong disagreement between the equa- 
tion of state obtained using the renormalized jellium 
model and the PBC theory. In the case of cylindrical 
polyions the disagreement is exacerbated by the fact that 
the ideal gas contribution to the equation of state, Eq. 
(|30p . vanishes in the limit of L — ► 00 considered in this 
work. For small 77, the behavior predicted by the renor- 
malized jellium model is more realistic than that of the 
PBC. 



V. CONCLUSION 

Starting from a mean-field description in which a dis- 
persion of N spherical or rod-like polyions is treated using 
a A-body Poisson-Boltzmann theory [29( , we have intro- 
duced the simplification of a homogeneous background 
to include the contribution of other colloids to the static 
field created by a tagged colloid. The charge of this back- 
ground is consistently renormalized to coincide with the 
effective charge governing the far-field potential. This re- 
sults in a non-trivial density dependence of the effective 
colloidal charge, which directly enters into the equation 
of state through a simple analytical expression. The good 
agreement observed between the pressure calculated us- 
ing the renormalized jcllium and the Monte Carlo simu- 
lations confirms the relevance of the renormalized jellium 
model for theoretical and experimental purposes and pro- 
vides an alternative to the Poisson-Boltzmann cell ap- 
proach. Furthermore, we note that the effective charge 
calculated using the renormalized jellium model should 
be more relevant for the study of the effective interac- 
tion between the colloids than its Poisson-Boltzmann cell 
(PBC) counterpart. This is particularly the case since at 
finite colloidal density, the DLVO potential arises natu- 
raly within the jellium formalism, while it has to be in- 
troduced extraneously within the PBC. In this work, we 
have left untouched the question whether the pair poten- 
tial calculated using jellium is a potential of mean-force 
or an effective pair potential (following the terminology 
of [![). Further work is required to answer this question. 

In a cylindrical geometry, the present approach implic- 
itly subsumes an alignment between infinite rods -which 
is also a prerequisite for the analysis of 0- but contrary 
to the crystalline structure underlying the introduction 
of the cell model, we consider here systems with no po- 
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sitional order for the rods. As for spherical colloids the 
pressure in both approaches is in good agreement up to 
relatively high densities, whereas the effective charges dif- 
fer significantly. We have also shown that the scenario 
for counterion condensation is similar to that of the cell 
picture. 

Our approach, which is best suited to describe systems 
with low macro-ions densities, may be easily extended 
to the case of asymmetric electrolytes. One interesting 
aspect of the renormalizcd jcllium is that the descrip- 
tion of colloidal mixtures (macro-ions with different sizes 
and charges) appears to be as straightforward as for the 
mono-disperse systems reported here. This is an impor- 
tant difference with the cell approach, which cannot be 
easily extended to such systems. 



Among the possible refinements, it is possible to con- 
sider an in-homogeneous jcllium with, again, rcnormal- 
ized charge. This should allow to extend the relevant 
range of densities where the model holds. Another in- 
teresting extension deals with the derivation of electro- 
kinetic properties. Work along these lines is in progress. 
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VII. CAPTIONS 

Fig. {l}: The effective charge as a function of back- 
ground charge for 77 = 10 -2 and -Zbare^s/a = 4 (spherical 
colloids, no salt). The physical solution Z c s — Zb ac k to 
the problem is the point of intersection with the first 
bisectrix (see inset, where a magnification of the relevant 
part of the main graph is shown). 

Fig. @: The effective charge (or equivalcntly, the 
background charge) as a function of bare charge for 
77 = 10~ 2 (spherical colloids, no salt). The line has slope 
unity to emphasize the initial "Debye-Hiickel" regime. 

Fig. ([3]): Comparison between the effective charges 
within the Poisson-Boltzmann cell (PBC) and the 
renormalized jcllium model, as a function of the volume 
fraction. Here Zbare^s/a ~ 6 (spherical colloids, no salt). 

Fig. (fJJ: Pressure as a function of volume fraction 
within the cell and the renormalized jellium model, on 
a log- log scale. Here, Z^Jb/cl = 6 (spherical colloids, 
no salt). The inset shows the same data on a linear scale. 

Fig. ([5]): Effective charge as a function of volume 
fraction within the PBC and the renormalized jellium 
model, for Abarc^B = 1- The inset corresponds to the 
saturation regime where Xbare — ► 00 (rod-like colloids, 
no salt). 



Fig. (JTUJ) : The effective charge for spheres as a 
function of volume fraction within the PBC and the 
renormalized jcllium model for n TCS a = 1 in the satu- 
ration regime Zbarc — ► 00. The inset shows the same 
quantity on a linear scale. 

Fig. (HJ): The osmotic pressure as a function of 
volume fraction within the PBC and the renormalized 
jellium model in the saturation regime for n Tes a = 1. 
The inset shows the same data on linear scale (spherical 
colloids). 

Fig. (|12p: Comparison of the PBC and the renormal- 
ized jellium equations of state with the one obtained 
in Ref. 32[ from the Monte Carlo simulations. Here, 
the macro-ion volume fraction is i] = 8.4 1(T 3 while 
Z haie £ B /a ~ 21.45. 



Fig. (fT5|) : The effective charge for highly charged 
cylindrical colloids (saturation regime) as a function of 
volume fraction within the PBC and the renormalized 
jellium model, for K rcs a = 1. The inset shows appearance 
of a minimum in the presence of salt. 

Fig. (|14|) : Osmotic pressure as a function of vol- 
ume fraction within the PBC and the renormalized 
jellium model for n rcs a = 1 and Xbare^B = 2. The inset 
shows the same data on linear scale (cylindrical colloids). 



Fig. (|6|): The effective charge as a function of the 
bare charge for different values of volume fractions 
(renormalized jellium model for charged rods). The 
inset shows a magnification of the main graph in the low 
charge regime. The present scenario is exactly that of 
the Manning-Oosawa countcrion condensation occurring 
in the cell model. 



Fig. ([7]): Manning radius Rm versus packing fraction 
for Abarc^s = 4.2. Extremely low densities have been 
considered to see the predicted power law dependence 
Rm oc r/" 1 / 4 , sec Eq. p^|) . The upper inset shows 
that the bare charge dependence of Rm also follows 
the form given by Eq. (jT5J) . The lower inset shows 
q(r)£B/e as a function of distance from the rod axis 
on a linear-log scale: as expected, the inflection point, 
indicated by the vertical arrow, coincides with the point 
where q{R.M)^B/e = 1. 

Fig. ©: Pressure as a function of volume fraction 
within the PBC and the renormalized jellium model, for 
both a moderately charged, and a highly charged rods 
(saturation limit), without added salt. The inset shows 
the same data on a linear scale. 



Fig. @: Ratio between k and K Tes as a function of 
volume fraction for n vcs a = 1 (saturation regime and 
spherical colloids). 
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